Birth and growth of cavitation bubbles within 
water under tension confined in a simple synthetic tree 
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Water under tension, as can be found in several systems including tree vessels, is met ast able. 
Cavitation can spontaneously occur, nucleating a bubble. We investigate the dynamics of spon- 
taneous or triggered cavitation inside water filled microcavities of a hydrogel. Results show that 
a stable bubble is created in only a microsecond timescale, after transient oscillations. Then, a 
diffusion driven expansion leads to filling of the cavity. Analysis reveals that the nucleation of a 
bubble releases a tension of several tens of MPa, and a simple model captures the different time 
scales of the expansion process. 



The dynamics of cavitation is a well studied process 
in free liquids (cavitation in the wake of propellers for 
instance), while cavitation in a microscopic confinement 
only recently received more attention. One important 
area of research is the transport of sap in trees and the 
potential blockage of this transport due to cavitation: 
these events can stop circulation of water in the vessels 
of real tree^i^, or in synthetic ones that are used for 
microfluidic flow transport driven by evaporation^—. 

The stability of water at negative pressure, i.e. under 
tension (see^ for a review), can be tested with hydrogels 
where the tension is fixed through the relative humidity 
of the surrounding mediurr^ii. However, such a method 
features only the final equilibrium state, and does not 
reveal about the cavitation process itself. The dynamics 
of cavitation involves fast time scales that are not known, 
apart from the fact that they are suspected to generate 
vibrations (noises in trees-). 

In this study we resolve the dynamics of cavitation 
of water inside hydrogel cavities. Cavitation can occur 
spontaneously or be triggered by a disturbance created by 
a laser pulse. The dynamics of the nucleated bubble and 
hydrogel cavity are recorded with a high-speed camera 
or a photodiode. 

We use pHEMA hydrogels of the same chemical com- 
position as in^''', but with a periodic array of holes or cav- 
ities, created using the novel experimental setup shown 
in Fig. la. An Nd:YAG laser (A = 532 nm, = 6 ns) 
is coupled to an inverted microscope. The laser pulses 
are focused by a 40x water immersion microscope objec- 
tive inside a thin layer of partially polymerized hydrogel 
that has been precured with UV light (dimensions ~ 1.5 
xl.5 cmxl mm). The hydrogel sample is bounded by 
two microscope cover slides separated by PDMS spacers. 
A single cavity is created after a laser-induced optical 
breakdown inside the hydrogel sample. The system is 
coupled to an automated translation stage that allows 
the design of a regular array of holes. The shape of the 
cavities strongly depends on the degree of polymerisation 
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FIG. 1: (a) Experimental set-up. (b) Microcavities in the 
hydrogel created by laser pulses with low UV pre-curing, (c) 
Microcavities when the hydrogel is pre-cured for a longer time. 



of the hydrogel: while they are smooth and spherical for 
low pre-polymerisation, they tend to form "cauliflower" 
shapes when more cured (Fig. [Il3,c), in a manner very 
similar to what is observed when forcing air in a poly- 
acrylamide gel^. Finally, the hydrogel is fully polymer- 
ized (Fig. lb). To fill the the cavities with water, the 
hydrogel is boiled with water for a few days, and soaked 
in degazed DI water. 

As the hydrogel samples are dried, water is set un- 
der tension because of evaporation. Cavitation sponta- 
neously occurs when tension is strong enough, or can be 
directly triggered with laser pulses, focused in the hydro- 
gel, in vicinity of water filled voids. The gel is sometimes 
slightly altered (apparition of grey spot. Fig. [2^). We 
speculate that a mild pressure wave is created, travels 
over the microcavity and thereby increases the liquid ten- 
sion for a brief moment. This mechanical disturbance in 
the hydrogel is able to trigger cavitation in neighboring 
microcavities, if under tension. The dynamics of cavita- 



tion is recorded with a high-speed camera at frame rates 
of up to 432,000 frames per second, with an exposure 
time of 1 /is. In addition, the central part of the image 
can be projected simultaneously onto a photodiode with 
a rise time of 0.03 /is and bandwidth of 12 MHz, pro- 
viding a measurement of the mean intensity in order to 
access faster dynamics. 

Fast time scales Our recordings reveal that a single 
bubble is nucleated, and that its volume quickly reaches 
a finite value, within 1 or 2 microseconds (Fig. [2^-c). 
We captured the first image of the nucleated bubble af- 
ter a delay of 1 /is after arrival of the laser pulse. If these 
images are motion blurred due to the exposure time of 
1/is, they show a transient non-spherical shape, see sec- 
ond image of Fig 2b. On the photo-diode signal we record 
transient high-frequency oscillations at typically several 
MHz, that we interpret as oscillations of the radius in Fig 
2c, with an initial velocity as fast as 30 m/s. We assumed 
the variation in light intensity to be proportional to the 
apparent bubble area on images. 

The following images of the high-speed recording show 
that the bubble after nucleation reaches a temporary 
equilibrium radius and then slowly grows. This is in 
contrast to cavitation nucleation in a free liquid, where 
the bubble expands to a maximum size and then rapidly 
collapses due to the pressure of the surrounding liquid. 
Experiments with an observation from the side also evi- 
dence that bubbles rise due to their buoyancy (Fig. [2]i), 
meaning they are not attached by a contact line and are 
quickly spherical. 

The instantaneous volume of the bubble (Vb) and of 
the cavity (Vc) are obtained by extracting an effective 
radius from image analysis (Fig. [3|), consisting in trac- 
ing iso- values lines on grey levels bitmaps. Interestingly, 
the volume of liquid Vi = Vc — H,o just after nucleation 
decreases. This is due to the compressibility of water, 
which has an adiabatic bulk modulus Ki = 2.2 GPa 
at 20° C. Therefore, the variation of liquid volume AVi 
(= AVc — Vb^o) gives information on the pressure change 
Ap in the liquid: 



Ap = Ki 



AVi 

Vi ' 



(1) 



We performed this analysis for a total of thirteen bub- 
bles, see Fig. [3l showing that the released pressure 
amounts to values ranging within the tens of MPa range 
(30 MPa ±16 MPa). Due to the measurement error in 
the image analysis, the estimated pressure presents a lot 
of uncertainty, but is comparable to that measured for 
cavitation in static hydrogels with the same composition 
(22 MPa)^'^, and to other experimental methods^. Note 
that the pressure released by triggered cavitation is of 
the same order than with spontaneous cavitation (Fig 
3d). Indeed, the laser perturbation is not very large, and 
in order to trigger cavitation we had to wait for the pres- 
sure to be close to the spontaneous threshold, i.e. when 
cavitation occurs spontaneously on nearby cavities. 
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FIG. 2: Image sequences showing the birth of a cavitation 
bubble (a) triggered by a laser pulse aimed away from micro- 
cavity, here on the left (b) spontaneous without triggering, 
here under 40x objective, (c) Photo-diode recording on a cav- 
ity with a radius of 40 /im: Estimated apparent bubble radius 
as a function of time, normalized by the bubble radius just 
after nucleation, for a single experiment. Red : experimental 
data ; black: slightly filtered signal. Insert: recorded relative 
variations of the photodiode signal from which the radius is 
estimated, (d) Long term rise of a bubble in a vertical sample. 



From this pressure, the stored elastic energy of wa- 
ter can be evaluated to be Ep = —^ApVb. If all 
this energy was released into kinetic energy, writing 

= ^piVbRb for a bubble in liquid^^ of density pi, we 
could expect bubble oscillations with maximum velocity 
Rb = {—Ap/SpiY^'^ ^ m/s, and a characteristic time 
r = Rb/Rb of 0.1 /is only (for a cavity such that i?5=10 
/im, i?c=40 pm). This time is of the same order than 
the observed initial relaxation time of the light intensity 
signal (0.3 ps). The system shows oscillations that are 
rapidly damped, probably because of dissipation in the 
gel surrounding the microcavities, besides viscous dissi- 
pation in water or acoustic emissions. 

One may wonder about the content of the cavitation 
bubble. After cavitation, the bubble quickly fills with 
water vapor, which has a diffusion coefficient I^vap = 
2.4 X 10-^ m^s-^ With Rq the initial size of the bubble, 
the typical vapor filling time is Rq/D^^p ^ 10~^ s, thus 
longer than initial dynamics. Then it also probably fills 
with air, that was dissolved in water and in the hydrogel. 
These processes might slightly change the pressure inside 
the bubble, but cannot account for the continuous post- 
nucleation expansion of the cavitation bubble. 

Long time scales. Indeed, after its fast expansion to 
the volume Vb^o^ the bubble is not static but slowly grows. 
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FIG. 3: Evolution of (a) the bubble radius Rb and (b) the cav- 
ity radius Rc. (c) We infer that the bubble volume increase 
is larger than cavity volume increase, meaning that the liq- 
uid volume decreases, (d) Estimation of the tension that is 
relaxed, for triggered (full circles) and spontaneous bubbles 
(open circles). 



until it completely fills the cavity, after a time of approx- 
imately 100 s. To quantify this growth, we analysed the 
evolution of the volumes of the cavity and of the bubble 
versus time (Fig. d]). The bubble expansion Vb — Vb^o ver- 
sus time follows a power law over more than 5 decades of 
time, with exponents in the range 0.35 to 0.45 (Fig. \^). 

Neglecting evaporation, there are only two ways to 
make the bubble grow: increase the volume of the cavity 
or transfer liquid from the cavity to the gel. The analysis 
of Vc{t) shows that the size of the cavity indeed increases 
with time (Fig. [Sb). However, since the bubble eventu- 
ally fills the whole cavity, an outwards flow of water must 
occur. 

We developed a simple model to account for the power- 
law expansion of the cavitation bubble. In the following, 
we will only describe the growth of the bubble by flow in 
the gel and hypothesize a constant cavity volume. 

When cavitation occurs, the pressure in the liquid 
suddenly relaxes from a homogeneous negative value 
pi = — Ap, to a local value pi = Pg — 2a /Rq with 
2a /Rq the capillary pressure due to water surface ten- 
sion, and Pg the gas pressure inside the bubble. We 
then have \pi\ <C Ap. This locally changes the osmotic 
pressure H = P — pi in the gel near the surface of the 
cavity, where P is the external isotropic stress (here be- 
ing constant, equal to atmospheric pressure). Incidently, 
a gradient of osmotic pressure VII develops across the 
gel, entailing water fluxes q which follow Darcy's law 

^ ~ ~^^^^Pi ~ ^77^ ^^^^ ^ permeability, /i the 
viscosity, (j) the polymer volume fraction. 

The conservation of mass implie^^: ^ = 0V • q. A 
change in volume fraction results in a change in osmotic 
pressure, and the last equation can be rewritten as: 



^ = i^sw V. 
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FIG. 4: (a) Long time growth of the bubble volume (red line: 
photodiode signal after oscillations, blue circles: camera) and 
of the cavity (14: dashes, AVc open squares). The bubble 
volume Vh is fitted with Eq. [3] (continuous blue line), (b) 
Diagram showing the origin of the growth: diffusion of water 
in the hydrogel. (c) Growth exponent for the bubble size for 
triggered (full circles) and spontaneous bubbles (open circles) . 



where iCgw = is the swelling modulus of the gel, 

describing how the polymer network of the gel swells or 
shrinks in response to water pressure changes^^. Eq. [2] 
is similar to the one developed for the sedimentation of 
a colloidal solution^ ^. Assuming low variations of per- 
meability, we have approximated the right hand side by 
I^V^n. The water content and osmotic pressure there- 
fore follow a diffusion equation, with a diffusion coeffi- 
cient D = K^^k/ji. The value^^ of this coefficient is 
around D ~ 10~^^ m^s~^. 

The diffu sion o f pressure entail fluxes that follow q = 
k/ (ix Ap/ \fnDi^ over the cavity surface S = 47ri?^. The 
transferred volume is therefore dvi/dt = qS^ so that un- 
der the hypothesis of a constant cavity volume Vc^ the 
bubble volume departs at long times scales from the ini- 
tial volume Vq (given by Eq. [T]) along: 



(3) 



(2) 



with n= 1/2, and ta = (l^^j ' which represents 
the time for which the bubble volume has doubled. 

Fitting experimental curves of bubble growth of Fig. 
[3^ (for times t<ls), provides ta and exponent n values 
displayed on Fig. |4t. The value of n is close to the 
predicted value 1/2, except for the largest characteris- 
tic time ta- The observed range of t^ is also in excel- 
lent agreement with the expected value, predicted using 
the negative pressure measured in the first section: with 
Ap = 20 MPa, i^sw = 100 MPa, D = lO'^^ m^s'^ and 
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values extracted from one experiment of Rc = 23/im and 
Ro = 6.8/im, one finds ta = 0.08 s. The final filling of 
the bubble, when = Rc can be deduced from Eq. [3l 
and occurs at a time: 

where we have neglected Vb in front of Vc. Using the 
same numerical values as above, one finds tf = 115 s, in 
very good agreement with the experimental observations. 
Note that the cavity volume is not constant in the exper- 
iment, so that the diffusion model should describe pri- 
marily the transferred liquid volume. Preliminary analy- 
sis of this volume provides larger growth exponents: this 
phenomenon is still under investigation. Also another 
challenge is to understand why, surprisingly, the curves 
of Fig. [4|i, show that there is no detectable transfer of 
fluid to the gel during the first 10 ms after cavitation, 
since Vb{t) and AVc{t) initially differ by a constant off- 
set. Then, an outward flow proceeds, which results in gas 
filling of the cavity. 

In conclusion, we report in this Letter the first direct 
observations of spontaneous and triggered cavitation of 
water under static tension in confined geometries. The 
nucleated bubble initially expands at a sub-microsecond 
timescale, driven by a tension of tens of MPa. These 
extreme tensions result in visible water volume changes. 
These extereme tensions lead to uncommon facts in mi- 
crofluidics: very fast velocities (30 m/s) and visible liquid 
volume changes associated with inertia, as can be seen on 
the initial damped oscillations of the bubble. Work is un- 
der way to fully characterize these oscillations. Also, the 
cavitation events studied here show a milder dynamic as 
compared to a free liquid with successive collapses. Thus, 
cavitation in water under tension leads to nucleation of 
bubbles but not to the destructive effects observed during 
cavitation bubble collapse. The presence of hydrogel in 



trees has been demonstrated recently^^ and could thus 
have a role in the prevention of damage by cavitation 
events. 

We also investigated the following slow growth of the 
bubble, which is limited by liquid flow in the surrounding 
gel, rehydrating it. Besides being useful for the design of 
hydrogel-made microfluidic devices, this fact leads to two 
remarks. 

First, diffusion in a hydrogel is a complex process: the 
coefficient D arises from a subtle balance between the 
elasticity and the permeability of the gel (Eq. [2]), which 
change when it dries or hydrates. Since the hydrogel be- 
comes softer as it re wets, this might slow the diffusion 
process as the cavitation bubble expands, explaining the 
fact that the growth exponents seen on Fig. HJare always 
lower than expected from normal diffusion. Another chal- 
lenge is to explain the origin of the latency time of about 
10 ms after cavitation when no flow towards the gel oc- 
curs. 

Second, the release of tension propagates in the gel 
when is re wetted by water, and this effect should be felt 
in closeby cavities. Indeed, when focusing the laser in the 
middle of three cavities at the same distance from each 
one, we noticed that the last triggered cavitation required 
much more laser energy than the first one. This inter- 
action might stabilize neighboring cavities, but only for 
limited times, since diffusion eventually smoothes out the 
temporary pressure peak. Such a phenomenon suggest 
that cavitation propagation can be hindered by negative 
correlations. Applied to the context of sap transport, 
we hypothesize that this phenomenon may play a role in 
trees to regulate tension in the xylem. 
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At 20° C, The adiabatic and isothermal compressibilities of 
water are very close and differ less than 1 percent. Their 
variation with pressure is also small: about 2.5 percents 
each 10 MPa. 

Conservation of mass for the polymer implies (l)SV = est 
for an elementary hydrogel volume 6V. Its variation is 
thus: ^^7- — while conservation of mass for the 

at (p at 

fluid implies: ^ = ^ = -^V^V • q, which gives the 
formula in the text. 



This swelling modulus under constant isotropic stress can 
also be expressed as K^^ = V^^ \ being the total hydro- 
gel volume. It is equivalent to a compression bulk modulus, 
but with opposite sign, as a polymer network inflates when 
increasing the solvent pressure. 

This value is deduced from measurements of the relax- 
ation time of a drying hydrogel. For our pHEMA hydro- 
gel, we estimate K^^ 100 MPa from Flory's theoryi^, 
so the estimated hydrogel permeability is very low: kf/n = 
10"^^mVPa.s. For comparison, with Polyacrilamide hy- 
drogel: Tanaka et al^^ find L> - 4 x 10~^^ m^s"\ This 
diffusion coefficient is similar to the present pHEMA hy- 
drogel, but originates from very different values of elas- 
ticity and permeability, since they have k/fi = 1/f = 
5 X 10~"'^^m^/Pa.s, but a much softer gel. 



